Pranck-Condon Factors as Spectral Probes of Polaron Structure 



OS 
Os 
Os 

G 

G 






c3 



i 

-a 
g 
o 
o 



> 

oo 
m 

CN 

O 

as 

Os 

^— > 



I 

G 

o 
o 



X 



David W. Brown 1 , Aldo H. Romero 2 , and Katja Lindenberg 1,3 
Institute for Nonlinear Science, 
University of California, San Diego, La Jolla, CA 92093-0402 

2 Max-Planck Institut fur Festkorperforschung, 
Heisenbergstr. 1, 70569 Stuttgart, Germany 

3 Department of Chemistry and Biochemistry, 
University of California, San Diego, La Jolla, CA 92093-0340 

(February 7, 2008) 



We apply the Merrifield variational method to the Holstein 
molecular crystal model in D dimensions to compute non- 
adiabatic polaron band energies and Franck-Condon factors at 
general crystal momenta. We analyze these observable prop- 
erties to extract characteristic features related to polaron self- 
trapping and potential experimental signatures. These results 
are combined with others obtained by the Global-Local vari- 
ational method in ID to construct a polaron phase diagram 
encompassing all degrees of adiabaticity and all electron- 
phonon coupling strengths. The polaron phase diagram so 
constructed includes disjoint regimes occupied by small po- 
larons, large polarons, and a newly-defined class of compact 
polarons, all mutually separated by an intermediate regime 
occupied by transitional structures. 

PACS numbers: 71.38. +i, 71.15.-m, 71.35.Aa, 72.90. +y 



I. INTRODUCTION 

The theory of polarons has undergone an evolution in 
recent years that has substantially improved our ability 
to put solid, quantitative accuracy to matters that have 
heretofore enjoyed only semi-quantitative estimation or 
qualitative characterization. This can be said in view 
of a convergence of results 0] that has been found in a 
number of independent and high-quality methods that 
have been brought to bear in particular on the analy- 
sis of the Holstein molecular crystal model [|p]. Im- 
portant among these methods are cluster diagonalization 
[0-M , density matrix renormalization group B , quantum 
Monte Carlo p-(14[, and certain variational approaches 
[0,15 -25 1 . Though quite distinct in their conception and 
implementation, these methods have all been found to 
be in deep and broad quantitative agreement over wide 
regions of the polaron parameter space. 

Our own effort in this area has relied mainly upon the 
Global-Local variational method, certain results of which 
will figure in the present work. A significant part of this 
effort has dealt with the problem of developing a reliable 
and interpretable polaron phase diagram on which can be 
clearly delineated the distinct regions of the polaron pa- 
rameter space where distinct classes of polaron structure 
may be found. In the course of this development, some 



familiar notions that have become part of the polaron 
lore have had to be revised, including the real physical 
nature of the large polaron p3,E3[,E6[ and the meaning of 
self-trapping in D dimensions p^ , |28| | . 

Here, we continue to be concerned with the polaron 
phase diagram, but in a manner and regime that are 
complementary to what has been already been devel- 
oped. For practical and formal reasons, the utility of 
the Global-Local variational method deteriorates signif- 
icantly when the fundamental electron transfer integrals 
are small and the electron-phonon coupling is weak; as a 
practical matter, this limitation excludes a sufficient por- 
tion of the non-adiabatic regime to preclude a meaningful 
assessment of the self-trapping transition there. Neces- 
sarily, therefore, what we have been able to say about 
the polaron phase diagram in the non-adiabatic regime 
has been based on extrapolations from more adiabatic 
behaviors. 

The non-adiabatic regime is important to many 
narrow-band systems and particularly to molecular crys- 
tals for which the Holstein model was originally formu- 
lated p|,p|pS| -jjjjj]. Polaron properties in the non-adiabatic 
regime generally depend quite smoothly on the basic sys- 
tem parameters, without the relative abruptness that 
tends to emerge in the adiabatic regime, and the low 
orders of perturbation theory, either weak-coupling or 
strong-coupling, tend to do a reasonable job of capturing 
most behaviors. Paradoxically, perhaps, it is this relative 
unremarkableness of the non-adiabatic regime that raises 
some of the questions motivating our study; in particu- 
lar, how the dramatic character of the self-trapping tran- 
sition that is so obvious at high adiabaticity dissembles 
into relative obscurity, and how, as a practical matter, 
its lingering presence may be recognized in observable 
polaron properties. 

We approach this problem through the use of the Mer- 
rifield variational method jl8| , |34| . The Merrifield method 
can be viewed as an antecedent to the Global-Local 
method in that it is the first in a sequence of increas- 
ingly refined variational methods leading to the Global- 
Local method. Although the Merrifield method suffers 
some very characteristic limitations that restrict its use- 
fulness as a tool for implementing polaron theory at gen- 
eral points in the polaron parameter space, it is at its 
best in the non-adiabatic regime where computation by 
the more general Global-Local method becomes difficult, 



and is thus well-suited to the present task. Moreover, 
owing to its relative simplicity, it is possible to pursue 
results in D dimensions, and to obtain some degree of 
analytical guidance and support for numerical studies. 

In order to locate the self-trapping transition, we need 
to analyze observable physical properties that take on 
distinguishable asymptotic behaviors on each side of the 
transition and objectively locate a boundary that dis- 
criminates between these behaviors. Here, we focus on 
two properties that are particularly important to spectral 
studies, the polaron ground state energy (related in well- 
known ways to Stokes shifts) and Franck-Condon factors 
(related in well-known ways to oscillator strengths). As a 
function of the electron-phonon coupling strength g, the 
ground state energy generally exhibits a "knee" between 
distinct weak- and strong-coupling trends that can be lo- 
cated and followed in parameter space to develop a self- 
trapping line. Franck-Condon factors generally exhibit a 
distorted Gaussian dependence on the coupling strength, 
allowing the central peak region (weak coupling) to be ob- 
jectively distinguished from the outer tail region (strong 
coupling). 

We use the Holstein Hamiltonian gg] on a D- 
dimensional Euclidean lattice 

D 

H = - } J } J Jial(aH+ei + afZ-rJ 

n i— 1 

+ fiUJ ^ b li b n - ghu X! °S a "(4 + & ") J W 

ft n 

in which at creates a single electronic excitation in the 
rigid-lattice Wannier state at site n, and ot creates a 
quantum of vibrational energy Tiw in the Einstein oscil- 
lator at site n. The Ji are the nearest-neighbor elec- 
tronic transfer integrals along the primitive crystal axes, 
and the li are unit vectors associated with the primitive 
translations. The above model encompasses all Bravais 
lattices, with the different lattice structures appearing 
only in the relative values of the hopping integrals Ji. 
For simplicity in the following, we use terms appropriate 
to orthorhombic lattices in which conventionally i = x, y, 
or z; however, all results hold for lattices of lower sym- 
metry with appropriate transcription of these labels to 
those of the primitive axes. 

We use the following Fourier conventions for ladder 
operators (c* = a^tf) and scalars: 
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N- D / 2 J2^ H ct, (2) 



>4 = N- D j2^' ri 4> >4 = E e ~^" A £- ( 3 ) 



It is convenient in the following to characterize tunnel- 
ing strength in D dimensions in part through a parameter 



J — J2t Ji', when restricting discussion to isotropic tun- 
neling, we use the notation J = Ji, such that J = DJ in 
those cases. 

For the most part in this paper, we limit our discus- 
sion to the non-adiabatic regime, defined by the condition 
J /hco < 1/4. Polarons at such small J /too are quite nar- 
row since we know that the largest polaron in any dimen- 
sion (as characterized by the size of the phonon cloud) 
has a width of y / 2Ji/hw along the i axis J26|; since no 
Ji/hco is greater than 1/4 in the non-adiabatic regime, no 
polaron in this regime has a width greater than a lattice 
constant, even at vanishing coupling. Thus, the varia- 
tional space in which the problem is solved numerically 
need not be large in order to contain the complete po- 
laron. This ability to contain the present problem in a 
small real-space volume facilitates computation consid- 
erably. 

The Merrifield trial state may be written 

|*(k)) = N-°/ 2 J2 e "'™4 

ft 

xexp[-A^/ 2 E(Afe- l «-"&t-/i.c.)]|0) , (4) 
q 
(tt(*)| * (£')) = fo, , (5) 

in which the {A5} are the variational parameters speci- 
fying the coherent state amplitude in the phonon mode 
q. Though these trial states are delocalized and satisfy 
the appropriate Bloch symmetry condition, and thus any 
property measured in the "lab" frame is uniform over the 
lattice, the internal structure of these delocalized states 
is determined by exciton-phonon correlations that are es- 
sentially local in character. Here, that local character is 
such that the electronic component located at site n is 
associated with a "phonon cloud" centered on that site, 
determined by the set of lattice amplitudes {A£}. 

We evaluate the variational energy band as the expec- 
tation value of the Holstein Hamiltonian 

D 

E n = {^{R)\H\^{K)) =-J2 Me %Kz S+ t + e- iKi Sli) 

+N- D ha J J2\^\ 2 - N- D ghu J Y / (^+ % (6) 

wherein S^ is the Debye- Waller factor 

Sl = cxp[iV-^|Af| 2 (e ±l<? *-l)]- (7) 



Minimization of E K with respect to X~* leads to the 
self-consistency equations for X'i: 



gfiuj 
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Sf = e X p[7V-^|Af| 2 (co Sgi -l)] 



(8) 



(9) 
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(10) 



(11) 



where S*f and ±$f are the magnitudes and phases of 
the complex Debye- Waller factors S^. This shows the 



optimal A£ to be real, and establishes the "sum rule" 



7 , *rt — -V=o — 5 



(12) 



valid at any k and in any number of dimensions. 

When any particular Ji/hui — > 0, A| loses any depen- 
dence on qi and m, becoming "flat" in those variables. 
The Debye- Waller factor (Sf ) and phase ($f ) associated 
with that direction drop out of the problem and the real 
space phonon amplitudes X'i become completely localized 
along the i axis. Although a disparity among the relative 
magnitudes of several Jj can result in a polaron that is 
in respects "small" in certain directions and "large" in 
others, there is not a distinct self-trapping transition as- 
sociated with each Ji [E8| . This can be seen in the present 
context in the fact that the Debye- Waller factor S z asso- 
ciated with a vanishing J z does not approach e~ 9 , which 
would be expected of ID small polarons along the z axis, 
but a quantity that is characteristic of the 2D polaron 
structure in the two surviving dimensions, whether this 
be large-polaron-like or small-polaron-likc. 

The sum rule continues to be satisfied as dimensions 
are turned off or on (e.g., J z — > in three dimensions), 
since 



lim 



\ \^x t^y i^z ) 



im > X) 

y ^n L 4 {n x ,n y ,n z ) 



ria-.rbv .n. 



E a 



(K x ,Ky)r _ ST^ \( K 



) =9 

y ) 



(13) 



Thus, there is no need for dimension-specific formula- 
tion if dimensions are controlled through the tuning of 
{Ji/hw}. 

Owing to the symmetries ([ll]) and the periodicity of 
the reciprocal lattice, «3>f vanishes at the Brillouin zone 
center and everywhere on the Brillouin zone boundary. 
Consequently, we have certain special values that play a 
significant role in the following. Denoting the reciprocal 
lattice origin by and any of the most extreme points of 
the Brillouin zone by 7? (ft, = ±7r along each axis), we 
find 
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(14) 



(15) 



The zone-center phonon amplitudes (O) are well- 
behaved under all circumstances because the denomina- 
tors, similar to those of weak-coupling perturbation the- 
ory, are sums of bounded positive terms. The zone-edge 
relation (119), on the other hand, suggests the possibility 
of encountering large or divergent phonon amplitudes for 
phonon wave vectors near the Brillouin zone boundary 
if tunneling is sufficiently strong (J'/hui > 1/4) and the 
Debye- Waller factors {Sf} are sufficiently near unity (as 
generally occurs when electron-phonon coupling is suffi- 
ciently weak; see the ID examples below). This potential 
divergence is both a real physical phenomenon and a gen- 
erator of artifacts in the Merrificld method. 

The reality of the phenomenon is due to the resonance 
or near-resonance that can occur between the states of 
the one-phonon continuum and zone-edge states of both 
the free-electron and the self-consistent polaron when the 
energy gap between latter and the one-phonon contin- 
uum is small. This circumstance occurs in any number 
of dimensions when J /Tiui > 1/4 and g is small. Un- 
der these circumstances, only a very small amount of 
electron-phonon coupling is needed to produce intense 
interactions that flatten the outer polaron energy band 
(level repulsion) and create heavy phonon clouds strongly 
modulated by the character of the zone-edge phonons. 

The Merrificld method accommodates the nearness of 
the one-phonon continuum by producing strong distor- 
tions of the variational amplitudes of a qualitatively ap- 
propriate nature; the phonon amplitudes become highly 
focussed around a single phonon wave vector, in clear ap- 
proximation to the single-phonon quantum that consti- 
tutes the exact g — ► state. However, because the Mer- 
rifield state is not well-equipped to emulate the highly 
quantum mechanical character of such states, the energy 
balance central to the variational principle is distorted 
and the variational energy bound is raised. Consequently, 
rather than experiencing the expected strong repulsion 
from the one-phonon continuum, the outer-zone portions 
of Merrifield energy bands flatten relatively weakly and 
cannot be taken as appropriately indicative of the po- 
laron structure when J fhio > 1/4. 

Although one cannot rely upon the numerical values 
of Merrifield band energies influenced by resonances with 
the one-phonon continuum, it is nonetheless true that the 
general character of the variational lattice state responds 
to such resonance effects in a reasonably appropriate way, 
provided that J /Tito < 1/4. This suggests that changes 
in the variationally-dctermined Franck-Condon factors, 
as very direct figures of merit for this general character, 
may reasonably indicate where the essential changes in 
polaron structure occur. It is thus that in the following 
we rely upon Franck-Condon factors as our primary di- 
agnostic of outer-zone polaron structure. Although this 
proves to be a very practical election, it is a choice that 
is in respects forced upon us by the limitations of the 
Merrifield method. As a choice that in the larger picture 
should be sufficient, but not necessary, other theoretical 
methods not so limited should find similar behavior in 



the band energies and other polaron properties near the 
band edge. 



II. THE ID CASE 

In principle, the set of equations (ph - (M) can be closed 
in the Sf and <Ff alone, greatly reducing the size of the 
self-consistency problem to be solved. This is of practi- 
cal advantage only in one dimension, however, since the 
reduction to quadratures involved in higher dimensions 
does not significantly facilitate computation. 

Replacing the summations in Eqs. (Q) and (|10|) with ID 
integrations, one arrives at the self-consistency equations 
first obtained by Merrifield: 



S" = exp(- 5 2 A K ) , 



(16) 



second order in g, and strong-coupling perturbation the- 
ory through first order in J/Tiuj. The latter is actually 
a shortcoming of the Merrifield method since important 
contributions from second order quickly dominate the 
first order of strong-coupling perturbation theory; how- 
ever, the first order is sufficient to properly determine 
that in the J/Hlu — > limit the "knee" in the dependence 
of E K on g at any k is given by 



d 3 E K 



3 J „ 

at — =0 



(20) 



The same differential criterion can be applied at finite 
J/tiu), which we will use to develop the phase diagram in 
Section IV below. 



III. THE FRANCK-CONDON FACTOR 



and 

A K : 



$ K = -g 2 A K (2JS K /hLo)sm($ K - - k) , 



fko 
y/[hto + 2JS' K cos($ K - k)] 2 - (2JS K ) 2 



(17) 



(18) 



Using ([H]) - ( |l8| ) in (g) yields the full set of variational 
phonon amplitudes such as shown in Figure [lj. 




FIG. 1. Sample surface showing the (real) variational am- 
plitudes \q in the ID case for J/Juu = 0.2 and g — 1. 

One may further obtain the energy-momentum rela- 
tion: 



E K =g 2 huj(A K -2[A K }3) 

-2JS K (1 - g 2 A K ) cos(k - $ K ) 



(19) 



The price paid for the compactness of this expression 
is the self-consistency condition that makes ( |l9| ) awk- 
ward to analyze; however, it can be shown that J19] ) 
agrees with weak-coupling perturbation theory through 



A quantity intimately related to the Debye- Waller fac- 
tors appearing in the self-consistency equations is the 
Franck-Condon factor 



F(,?) = |<^)|4|0)| : 



(21) 



This is one of many Franck-Condon factors associated 
with various transitions between correlated electron- 
phonon states. This particular factor characterizes the 
oscillator strength of the zero-phonon line associated with 
a transition between a free electron state of crystal mo- 
mentum k and the polaron it forms at the same k; these 
are direct transitions, resolved by crystal momentum. 
Though one is often concerned primarily with transitions 
near the Brillouin zone center, we will use the full k de- 
pendence of F(k) across the Brillouin zone, and particu- 
lar at the zone center and the zone boundary. 

Our principal interest in zero-phonon lines in the 
present context is in the possibility that they may offer 
an observable means for mapping the essential polaron 
features on the polaron phase diagram. Owing to the 
strong similarity between the Franck-Condon factor and 
the Debye- Waller factors intimately connected with the 
polaron effective mass, it is reasonable to expect that an 
analysis of the dependence of the Franck-Condon factor 
on model parameters should be able to yield the loca- 
tion of the self-trapping line. Also owing to the fact that 
the Merrifield method is at its best in the non-adiabatic 
regime and weak coupling, it is reasonable to hope that 
such an analysis would complement others made by other 
methods generally more accurate (e.g., the Global-Local 
method), but which deteriorate in quality at small J /Tvu) 
and g. 

The quenching of the zero-phonon line is an experimen- 
tal signature that has long been associated with the self- 
trapping transition. As a function of model parameters, 
this quenching occurs continuously as electron-phonon 
coupling is tuned from the weak-coupling to the strong- 
coupling regimes. There is thus some inherent ambiguity 



in the assignment of the point we associate with the self- 
trapping transition; our criterion here, as elsewhere, is to 
identify the self-trapping transition as the point of most 
rapid change in a property that takes on characteristi- 
cally different behaviors in the weak and strong-coupling 
regimes. In the case of the Franck-Condon factor, this 
criterion takes the form of an inflection point in the de- 
pendence of F(k) on g at fixed {Ji}. 

In terms specific to our variational development, 



FW = exp[-N- D J2\>4\ 



(22) 



which is just the exponential of the average number of 
phonons per mode in the phonon cloud. 
In the ID case, we find 

F(k) = exp {-.9 2 A K [1 - (2JS K /huj) cos(k - $ K )]} , 

(23) 

as shown in Figure || in selected cases. 
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(See Figure |[) It is clear from these that the Franck- 
Condon factor is independent of K at J /Two = 0. It is also 
clear that increasing J/huj from causes F(0) to broaden 
out to stronger coupling and causes F(tt) to narrow to- 
ward weaker coupling. 
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FIG. 3. Franck-Condon factors in D dimensions at the 
zone center (upper curves) and at the most extreme point of 
the Brillouin zone (m = ir) for DJ/Tiuj — — 0.25. 

Perhaps the most interesting behavior revealed in 
© - @ is that of F(tt) at J/Tiuj = 1/4. Using the 
fact that at this particular J /Tito value 

S n = exp[-g 2 (l - S n )- 3/2 } , (26) 

one can show that the leading dependence of S n on g is 



S" 1 ' ~ I 



9 



4/5 



(27) 



This in turn implies that in the same approximation, 



FIG. 2. Franck-Condon factors across the Brillouin zone 

in ID. J/Tiuj = 0.2, g = 0.25 - 1.5. 

It is simple to show that in the absence of tunnel- 
ing F(k) = e~ 9 for all k in any number of dimensions. 
Thus, in the J /Two — > limit the self-trapping features 
associated with the Franck-Condon factor are found at 
g = l/v2, for any k. With increasing tunneling, this de- 
generacy is broken and these features fan out; the manner 
in which this occurs can be seen most clearly in the ID 
case, where 



F (k = 0) = exp 
F(n — 71") = exp 



1 - 2JS°/huj 
(1 + 4JS /Tico) 3 / 2 
2 1 + 2JS n /huj 
(1 - AJS^/hiof/ 2 



(24) 
(25) 



F(n) ~ exp[- 



,4/5i 



(28) 



(See Figure [|.) This singular behavior in the Franck- 
Condon factor at J /Two — 1/4 suggests that the self- 
trapping feature identified by an inflection point in g 
moves to g = at J /Two = 1/4. Indeed, on can show 
that for J '/Two = 1/4 — e, the Franck-Condon factor re- 
tains an initial finite negative curvature in g, suggesting 
that there exists a proper inflection point at finite g. 
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FIG. 4. log(— In F (/?)) vs. log g, allowing the power of g in 
the exponent of the Franck- Condon factor to be ascertained. 
Curves are coded to correspond to Figure 0; zone center (lower 
curves) and zone edge (upper curves) in D dimensions for 
DJ/Huj = 0, 0.05, 0.1, 0.15, 0.2, and 0.25. 

Thus, all these considerations suggest that the self- 
trapping feature in F(0) should should shift more-or-less 
steadily from g = l/v2 at J/fiuj — to stronger cou- 
pling with increasing J/ftw, while the self-trapping fea- 
ture in F(n) appears to shift from the same g = l/v2 to 
weaker coupling with increasing J/Hlo, arriving at g = 
at J/huj = 1/4. 

The particular results above have been obtained for the 
ID case, which enjoys sufficient tractibility to admit some 
reasonably straightforward formal analysis. In higher di- 
mensions, numerical solution and analysis is generally 
more practical, though certain exceptions warrant special 
attention. Throughout the foregoing we have highlighted 
the particular crystal momentum values associated with 
the Brillouin zone center (re = 0) and the most remote 
corners of the Brillouin zone where all the crystal mo- 
mentum components take their maximum values (re = ff). 
We demonstrate in Appendix A that for these particu- 
lar re values in the isotropic cases in D dimensions, the 
dependence of the Franck-Condon factors on the dimen- 
sionality and tunneling strength is reduced to the single 
scaled variable D J. This implies that the particular re- 
sults shown in Figure [| hold not only in ID, but in 2D 
and 3D as well when parameters are scaled appropriately. 

Results can be obtained numerically for any degree of 
anisotropy and general re; however, the qualitative na- 
ture of the dependence on anisotropy can be inferred 
from Figure y without detailed analysis. Consider, for 
example, re = and J/fwo = 0.05: The results for 
quasi-2D scenarios with J x /h~ui — Jy/hcu = 0.05 
and < J z /fLui < 0.05 are contained between the 



DJ/huj = 0.15 and DJ/hui = 0.1 cases shown in 
Figure |[ Similarly, the results for quasi-lD scenar- 
ios with J x /huj = 0.05, < Jy/huj < 0.1 and 
Jz/Tiui = are contained between the DJ/hui = 0.05 
and DJ/hcj = 0.1 cases shown in Figure LI 

The tunneling parameters and the effective dimension- 
ality determined by them are not generally subject to any 
practical degree of experimental control; however, even 
greater changes in the Franck-Condon factor can be in- 
duced by changing the magnitude and/or the orientation 
of the crystal momentum re* probed. To the extent that it 
is possible to achieve some selectivity in the re's sampled 
in a particular experiment, it should be possible to induce 
controlled variations in the oscillator strengths associated 
with these Franck-Condon factors by, for example, vary- 
ing the orientation of the sample. Such predictable "wob- 
bles" in the intensities of appropriately-selective spectral 
probes constitute signatures of polaron structure that ex- 
ist only if both electron-phonon coupling and electron tun- 
neling are sufficiently great, without reducing the quasi- 
particle to the status of a "mere" small polaron. 



IV. PHASE DIAGRAM 

The overall character of the foregoing results can be 
summarized on a diagram of the polaron parameter space 
in which the loci of the knees in the polaron band energies 
and the inflection points of Franck-Condon factors play 
the role of rough phase boundaries separating distinct 
classes of polaron structure (see Figure H). These lines 
are accurately described by the simple relations 
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DJ/huj < 1/4 , (29) 

DJ/huo < 1/4 , (30) 

DJ/huj < 1/4 , (31) 

DJ/huj < 1/4 . (32) 



obtained in empirical fashion by noting the exact J = 
termini as discussed in prior sections and augmenting 
these with the simplest expressions in whole numbers 
that express the apparent trends in a quantitatively con- 
sistent way. The restrictions on (E9h and (Rlh are weak 
because such zone-center properties are well-behaved un- 
der the Merrifield method to J /fuv substantially greater 
than 1/4; however, the quantitative accuracy of the Mer- 
rifield method even at the zone center deteriorates signif- 
icantly with increasing J^/huj, warranting prudence be- 
yond the strictly non-adiabatic regime. On the other 
hand, the restrictions on ( p0| ) and ( |32| ) are strong be- 
cause zone-edge properties are strongly affected by the 
one-phonon continuum. 



These boundaries only roughly distinguish distinct po- 
laron regimes because the structural changes occurring 
in the non-adiabatic regime are quite smooth and broad, 
with changes in different aspects of polaron structure 
occurring with less synchronization than is seen in the 
adiabatic regime. This is seen clearly in the fact that 
the critical features of the band energies and the Franck- 
Condon factors are significantly separated in g in the 
J '/hio — > limit, and while trending similarly with in- 
creasing J'/Huj, remain well separated over the entire 
non-adiabatic regime. 
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FIG. 5. Polaron phase diagram in the non-adiabatic 
regime according to the Merrifield method in D dimensions. 
Diamonds: Location of inflection points in F(tt). Circles: Lo- 
cation of inflection points in F(0). Squares: Location of the 
knee in E°. Triangles: Location of the knee in E n . Solid line, 
Eq. B2. Long-dashed line, Eq. plj. Short-dashed line, Eq. pfl. 
Chain-dotted line, Eq. 



to be identified that appears to accurately describe the 
central trend of such clusters of data over essentially the 
entire adiabatic regime. 

Similarly, a zone-edge curve 



9n = 1 + V * 



hu 4 



-1/8 



(34) 



in ID has been identified that appears to accurately de- 
scribe the characteristic changes in the outer energy band 
that signal the onset of significant narrowing of the po- 
laron energy band, commencing the process that develops 
into the self-trapping transition with increasing coupling 
strength. 

It is telling to combine the empirical curves (|3^) and 
( p4| ) abstracted from our prior ID Global-Local analy- 
ses with the complementary curves (|2^), (pi]), and ( |3^ ) 
that follow from our present analysis by the Merrifield 
method. This comparison is presented in Figure @. 
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This disperse character of the collection of self- 
trapping-related loci is not an artifact of the Merrifield 
method, nor of the particular parameter regime nor of 
the particular physical quantities used to locate tran- 
sition effects. A similar and complementary dispersion 
has been found at somewhat larger J/fiuj in ID using 
the Global-Local method, based on the analysis of phys- 
ical quantities such as the ground state energy, kinetic 
energy, phonon energy, electron-phonon interaction en- 
ergy [p0| , effective mass plf , and electron-phonon corre- 
lation functions |E2). In such analyses, the self-trapping 
loci attributable to different zone-center physical quanti- 
ties have been found to cluster increasingly tightly with 
increasing J/fuv, permitting an empirical self-trapping 
curve 



(33) 



' / "' /=1 + ^ 



FIG. 6. Polaron phase diagram combining the empiri- 
cal curves (|29|), (|3l|), and (j32|) based on the present analysis 
employing the Merrifield method for DJ/Tilu < 1/4 with the 
complementary empirical curves H33) and (p4) based on inde- 
pendent analyses by the Global-Local variational method for 
J/huj > 1/4 in ID. 

It is clear that zone-center properties follow a common 
pattern of behavior, relatively indifferent to the crossover 
from the non-adiabatic into the adiabatic regime but for 
the possible role the latter may play in setting the scale 
over which the disperse self-trapping loci of the non- 
adiabatic regime begin to cluster more tightly toward 
the more sharply-defined trend in the adiabatic regime. 
This dispersity does not vanish suddenly at the crossover, 
but is apparent as well in the dispersity of similar loci 
found under the Global-Local method at small J jTiio still 
greater than 1/4. 



Further, it is interesting and surely no accident that 
the Merrifield zone-center lines (|2^) and (|3l]) intersect, 
that this intersection falls very nearly upon the Global- 
Local zone-center line (f53|), and that this cluster of inter- 
sections coincides well with the first appearance of dis- 
continuities in the solutions of the Merrifield method. It 
was, in fact, the set of such "critical points" (J c /huj, g c ) 
ascertained from our sequentially-refined variational cal- 
culations (Merrifield method |L8[ , Toyozawa method [19[ , 
Global-Local method |16[) that first led us to iden- 



tify the simple empirical curve (33) as a convenient 
method-independent approximant to the real physically- 
meaningful self-trapping line. The "critical" appearance 
of solutions near such points (J c /huj,g c ) is, of course, a 
methodological artifact, and we should view the intersec- 
tions of the several zone-center lines near J/huj w 0.9 as 
an artifact as well; the physical self-trapping line surely 
bends smoothly through this region, seamlessly joining 
the central trend of the non-adiabatic loci with the more 
sharply-defined trend that develops at higher adiabatic- 
ity. 

These zone-center results of the Merrifield method sug- 
gest an answer to one of the more empirical questions left 
open by our prior Global-Local analyses. Though quite 
accurate over a large range of J/huj, it has seemed un- 
likely that the dependence of the empirical curve gsr on 
y^J/hoj should continue unregularized all the way down 
to J/Hlu = 0. If we are to continue to regard the phys- 
ically meaningful gsr as representing a central trend in 
the inherently disperse set of self-trapping loci even as 
J/hui vanishes, the present results suggest that the lead- 
ing dependence of gsr on J/hco should not persist as a 
square root, but yield to a more pedestrian linear depen- 
dence 



gsr ~ 1 + a 



J_ 
ftio 



J_ 1 



(35) 



where a is a constant of order unity. 

It is likewise clear that the zone-edge loci follow a com- 
mon pattern of behavior, albeit one that is exquisitely 
sensitive to the crossover from the non-adiabatic into the 
adiabatic regime. Although the loci illustrated above and 
below this crossover are drawn from different physical 
properties (Franck-Condon factors from the Global-Local 
method are not available and the Merrifield band energies 
near the zone edge are not meaningful in and beyond the 
crossover regime) , they are closely related and reflect the 
same, if complementary, underlying physical behavior. 

We are led to view the results of the Merrifield method 
in the non-adiabatic regime and those of the Global-Local 
method in the adiabatic regime as mutually confirming, 
and describing one consistent set of self-trapping phe- 
nomena at all J/fruj. 

In the adiabatic regime, it is quite straightforward to 
view the two lines gsr and g^ as dividing the polaron pa- 
rameter space into a small polaron regime at strong cou- 
pling, a large-polaron regime at weak coupling, and an 
intermediate regime occupied by transitional structures. 



In the non-adiabatic regime, it likewise clear that there is 
an unambiguous small polaron regime at strong coupling; 
moreover, it is noteworthy that the non-adiabatic small 
polaron states are continuously deformable into the adi- 
abatic small polaron states without encountering tran- 
sition behavior in any basic polaron property, so that 
there is no formal distinction to be made between non- 
adiabatic and adiabatic small polarons. A complemen- 
tary observation can be made about large polarons; al- 
though large polarons as we have thus far characterized 
them reside in the adiabatic regime, there is no formal 
distinction to be drawn between large polarons at greater 
or lesser electron-phonon coupling strengths since these 
are continuously deformable into each other without en- 
countering any transition behavior. The intermediate 
regime can be defined in similar terms as that in which 
transition behavior is found in some basic polaron prop- 
erty at every point; for example, although the lines gsr 
(zone center) and <?at (zone edge) are only discrete curves, 
the domain between them is dense with similar curves 
associated with the occurrence of transition behavior at 
general k values. 

These observations lead us to consider the more darkly- 
shaded region of Figure @ at the weak-coupling end of the 
non-adiabatic regime. The transition line (32) associated 
with the zone-edge Franck-Condon factor and the more 
limited, qualitative information available in the Merri- 
field band-edge energies (|3C]) suggests that this regime is 
disconnected from both the small polaron regime and the 
large polaron regime in the sense of continuous deforma- 
bility as used above. Provided that this non-adiabatic 
weak-coupling regime is not dense with transition loci, 
which seems quite unlikely, it would appear that this 
regime is occupied by a polaron structure that is nei- 
ther "small" nor "large" nor of a transitional nature that 
would identify it with the intermediate regime. 

Indeed, some very basic polaron properties behave in 
qualitatively distinct ways in this insular regime, per- 
haps foremost being the polaron radius as given by the 
polaron Wannier function. Intuitively, one expects the 
polaron radius by any definition to decrease monotoni- 
cally with increasing electron-phonon coupling, and this 
is generally the case in the regimes we have here charac- 
terized as the small polaron regime and the large polaron 
regime. Even in the non-adiabatic weak-coupling regime 
now under discussion, the radius of the phonon cloud as- 
sociated with zone-center polarons has an initial width of 
order y / 2Ji/huj along the i axis Pq] , and decreases with 
increasing coupling strength. 

The polaron Wannier function, however, is a construct 
of the entire polaron energy band, being a superposition 
of polaron Bloch states of every k; as such, it is a localized 
state that can be viewed as energy band theory's own 
answer to the inverse problem of determining the identity 
of the localized quasiparticle whose dynamical properties 
are manifested in the polaron energy band. The real- 
space width of this state can be gauged in various ways, 
among them being a variance measure of the electron 



density within the polaron Wannier state. 

In present terms, we may construct polaron Wannier 
states from the trial Bloch states in the following fashion 

|$(n)) = --==$>-**■« !*(*)>' ( 36 ) 

from which we may construct the electron density profile 
as 

^=($(0)|4aH$(0)>- (37) 

Using this density, we may construct a variance tensor 

a % = E nr 3P f ( 38 ) 

in terms of which the spatial variance of the electron 
density in the polaron Wannier state may be given in an 
arbitrary direction. In the particular case of measure- 
ment along the x axis in three dimensions, for example, 
this result after summing over the y and z axes takes the 
form 
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With further manipulation it can be shown that 
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(41) 



The spatial variance of the localized electron density 
within the polaron is thus seen to be the average over 
all phonon modes and over polaron crystal momenta in 
the measurement direction of a mean square measure of 
the amount of distortion present in the the phonon am- 
plitudes along the measurement direction. 

Without further explicit calculation, this relationship 
provides a means of understanding what is distinct in 
the weak-coupling polaron behaviors found in the non- 
adiabatic and adiabatic regimes. In the adiabatic regime, 
the weak-coupling polaron band is nearly identical to that 
of the free electron in the inner Brillouin zone, but is 
strongly flattened in the outer Brillouin zone where the 
effects of interaction with the one-phonon continuum are 
severe. The phonon amplitudes exhibit strong changes 



in k which, through (41), are associated with broad po- 
laron Wannier states. With increasing electron-phonon 
coupling the severity of these K-dependent distortions de- 
creases, resulting in the narrowing of the polaron Wan- 
nier state. In qualitative terms, this narrowing trend is 
what is expected of large polarons. 

In the non-adiabatic regime however, quite a differ- 
ent situation is found. At weak-coupling, the polaron 
band is nearly identical to that of a free electron at all 
k, and in the limit of vanishing coupling is completely 
undistorted. The associated phonon amplitudes are not 
only very small, but are also very weakly distorted in k 
which, through (Ell) implies the complete localization of 
the polaron Wannier state on a single site as g — * 0. Con- 
versely, with increasing electron-phonon coupling, the 
presence of the higher-lying one-phonon continuum is felt 
more strongly at higher k, resulting in an enhancement of 
phonon amplitudes in the outer zone with which is associ- 
ated a slight flattening of the polaron energy band. This 
growth in K-dependent distortion results in a broadening 
of the polaron Wannier state with increasing electron- 
phonon coupling until a transition is made into the small 
polaron regime. 

Although such compact polarons n are straightfor- 
wardly understood, that they are completely localized 
in the weak-coupling limit and broaden with increasing 
electron-phonon coupling suggests that they be regarded 
as distinct from both the large polarons and small po- 
larons that dominate the outer reaches of the polaron 
parameter space. Moreover, that the compact polaron 
regime appears to be disconnected from both the large 
polaron regime and the small polaron regime by observ- 
able transition behavior suggests that such distinctions 
may be important to the clear classification of polaronic 
systems. 



V. CONCLUSION 

Our study of the Holstein model has focussed on the 
basic properties of observable zero-phonon lines in optical 
spectra; specifically, the polaron ground state energies 
that are principal determinants of the spectral position 
of such lines, and the Franck-Condon factors that are 
principal determinants of the oscillator strength of such 
lines. We have found in these results several properties 
that facilitate both the application of these findings to 



1 We use the term "compact" here to emphasize the com- 
plete spatial localization of polaron Wannier states that dis- 
tinguishes the weak-coupling limit of the non-adiabatic regime 
from that of the adiabatic regime. This usage is at least 
semantically consistent with that of the theory of nonlinear 
waves, where particle-like excitations with finite spatial sup- 
port have been termed "compactons" [B5H37|. 



experiments on real bulk materials and the interpretation 
of the experimental results in terms of underlying polaron 
structure. 

First, although quantitative results can be obtained for 
any polaron wave vector k, we have found that at the ex- 
tremes of isotropic polaron bands, at the Brillouin zone 
center and At its most remote corner, the dependence of 
the Franck-Condon factor on real-space dimensionality 
D and the elementary tunneling parameter J reduces to 
the single scaled parameter D J. This permits a straight- 
forward understanding to be had of how observations in 
systems of reduced effective dimensionality, for example, 
can be expected to be related to observations in bulk 
systems. 

Second, with some additional care, this general trend 
in the dependence of the Franck-Condon factors on di- 
mensionality can be seen to be similar to that which can 
be expected in the K-dependence of the Franck-Condon 
factor in a system of fixed dimensionality; for example, 
in changing the orientation of the probed wave vector 
from [1,1,1] to [1,1,0] to [1,0,0]. Such experimentally- 
controlled variations in the Franck-Condon factor (and 
its associated zero-phonon line) of a fixed system con- 
stitute a signature that can be associated with specific 
polaron structure. 

Third, beyond such quantitative characteristics we 
have found that the Franck-Condon factors at the Bril- 
louin zone center and at the extreme Brillouin zone 
edge constitute particularly direct means of revealing the 
changes in polaron structure that are associated with 
the self-trapping transition. That is, it appears possi- 
ble to map out the polaron phase diagram from surveys 
utilizing Franck-Condon factors alone. Through such 
considerations, here extended for the first time to the 
non-adiabatic weak-coupling regime, we have been able 
to complete a systematic appraisal of polaron structure 
spanning all regimes. The delineation of transition curves 
by means of the Franck-Condon factors has compelled us 
to distinguish a third kind of characteristic polaron struc- 
ture, the compact polaron, from the more familiar notions 
of the small polaron and large polaron. It is suggested 
that each of these three classes of polaron structure is 
separated from the others by an intermediate regime in 
which transition behavior is found in some basic polaron 
properties, but that this classification is essentially com- 
plete. 

Thus, suitably discriminating experimental studies of 
the detailed behavior of zero-phonon lines would appear 
to afford versatile, powerful, and interpretable means of 
deducing the structure of polarons in real materials. Suit- 
ably constructed surveys of such structure in a variety of 
materials should be capable of mapping such globally- 
important features as the polaron self-trapping line, pro- 
viding experimental tests of the proposed polaron phase 
diagram. 
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APPENDIX A: DIMENSIONAL SCALING 

In some of our illustrations of specific results we take 
advantage of a certain scaling property that holds under 
the Merrifield method at the Brillouin zone center and 
at selected points on the Brillouin zone boundary. The 
demonstration of this property utilizes the fact that the 
Debye- Waller phases $£ vanish at the Brillouin zone cen- 
ter and everywhere on the Brillouin zone boundary, and 
the fact that 
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as follows from M). 

Restricting discussion to the zone center (k = 0) and 
any of the most extreme corners of the Brillouin zone 
(k = 7?), we find that the fundamental variational ampli- 
tudes can be expressed in the form 
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Now further restricting to the isotropic case, we find 
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It is the reduction of the dimension- and J-dependence 
of the principal quantities to the simple combination D J 
that is responsible for the simplicity of our main results. 
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